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ABSTRACT 


The purpose of this thesis is to analvze and evaluate explosive cyclogenesis during 
the winter of the First Global GARP Experiment (Januarv to February 1979). Explosive 
cyclogenesis is defined as a decrease in the sea-level pressure at the rate of one mb per 
hour for a period of 12 h up to 24h. The European Centre for Medium-range Weather 
Forecasts (ECMWF) provided the revised analyses for evaluation and comparison of 
important cvclone properties in a sample of 13 explosive developing cases and eight 
nonexplosive cases. 

The specific parameters being examined include the static stability, low-level abso- 
lute vorticity, vorticity advection, eddy and mean modes of the vorticity transport. 
upper-level divergence, kinematic vertical velocities and the strength of the low-level 
baroclinity. These parameters are compared at the initial, 12 and 24 hour time periods 
as Well as the overall 24 hour average. The statistical relationships and magnitudes of 
these terms indicate the most significant physical mechanisms in explosive cyclone de- 
velopment compared to the nonexplosive storm group. 

The important outcome of these results is that the kinematic vertical velocity and 
the upper-level forcing mechanisms are statistically separable. The large values for the 
upper-level processes suggest that the upper tropospheric wave influence 1s most likely 


producing the stronger vertical motions. 
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IL. INTRODUCTION 


Rapid cvclone development in maritime regions of the world poses a particularly 
dangerous threat to today's mariner. These explosive cvclogenesis events are defined as 
low pressure svstems with a deepening rate of one mb per hour for a 12 h period or 
longer. They are primarily maritime events with a maximum occurrence in the western 
oceans. Current operational forecast methods have been unreliable in accurately pre- 
dicting the development of these types of storms (Sanders, 1987). The difficulty in 
forecasting these svstems is due to their often small scale, their rapid deepening rate and 
the lack of in-situ data in the ocean environment. 

Two recent examples of rapid cyclogenesis events discussed in the literature are: the 
Presidents’ Dav cvclone (18-20 February 1979) by Bosart (1981), Bosart and Lin (1984) 
and Uccellini et al. (1983, 1985); and the Queen Elizabeth II storm (9-11 September 
1978) by Gvakum (1983a,b). Anthes et al. (1983) and Uccellini (1986). 

Bosart and Lin's (1984) study of the Presidents’ Day snow storm illustrated the im- 
portance of the planetary boundarv laver processes and cold air damming to the east of 
the Appalachians for explosive cvclogenesis. Uccellim et al. (1984) discussed the im- 
portance of unbalanced flow in the subtropical jet near the eastern seaboard that re- 
sulted in coupled indirect vertical circulations. In this scenario, the low-level jet in the 
lower branch of this circulation aided in the development of precipitation even though 
the trough axis was farther upstream. 

The Queen Elizabeth II storm development was originally attributed to latent heat 
release by Gvakum (1983a). In this case, Gvakum asserted that the storm developed in 
an environment in which baroclinic support was confined to the lower troposphere. 
However, Uccellini (1986) pointed out that there were important upper-level baroclinic 
processes upstream of the rapidly developing cyclone. Twelve hours prior to the rapid 
cvclogenesis, a short wave trough was located 400 to 500 km. Associated with the 


trough was: 


e à polar jet marked by maximum winds of 65 m/s and significant vertical and hori- 
zontal wind shear; 


* positive vorticity advection and divergence at the 300 mb level; and 


e an intense front from the 300 mb level to the surface. 


The physical processes involved in the rapid deepening of cyclones are still being 
studied. There are two primary views on the physics of explosive cvclogenesis. The first 
approach by Bosart (1981) and Gyakum (19$3a.b) emphasized the low-level baroclinic 
environment and the convective bulk heating effects to explain the cyclone development. 
In addition, Rogers and Bosart (1986) documented that the explosive deepening cyclone 
resulted from a shallow baroclinic system with weak vertical stability near the center. 

The second approach by Uccellini et al. (1985) and Wash et al. (1988) concentrated 
on the upper tropospheric jet and wave structures. In this case, the rapid cyclone de- 
velopment was related to the interaction between diabatic and dynamical processes in 
the troposphere, including the strength of the upper-level jet and its ageostrophic circu- 
lation, and the static stabilitv of the marine boundary layer. 

A statistical study of many explosive and nonexplosive cyclones during the period, 
17 January to 23 February 1979 was conducted bv Smith (1986). This study of 14 ex- 
plosive and 14 nonexplosive cyclones concentrated on upper-level forcing mechanisms 
and boundarv laver characteristics during the initial rapid deepening periods. Diagnostic 
calculations for both storm groups were compared. His evaluation of storm environ- 
ment parameters suggested upper-level support was important in the deepening process. 
In the explosive cases, the positive vorticitv advection, eddv mode of vorticity transport 
and the upper-level divergence were significantly greater than in the nonexplosive cases. 
However. other important factors. such as the low-level baroclinicity, were not thor- 
oughly investigated. 

The intent of this research is to evaluate explosive maritime cyclogenesis using the 
new, revised First Global GARP Experiment (FGGE) analvses (Paegle, 1986). The main 


Objectives ol thiis thesis are to: 


e update the collection of explosive and nonexplosive cyclones of Smith (1986), using 
the enhanced FGGE analyses; 


e evaluate a number of cyclone diagnostic terms to study storm intensity, low-level 
baroclinicity, static stability and upper tropospheric processes for the explosive and 
nonexplosive cases: and 


e statistically determine dynamical differences between explosive and nonexplosive 
cyclogenesis. 

In this study, the quasi-Lagrangian or storm following budgets discussed by 

Johnson and Downey (1975a,b) are used to determine the cvclone's properties and 

transports. The cyclone development is related to the dynamical and thermodynamic 


processes in the upper and lower troposphere. The horizontal and vertical mass circu- 


to 


lation is directly proportional to the amount of low-level convergence and upper-level 
divergence. The mass budget analysis describes this vertical and horizontal transport. 
The important features of upper-level forcing that are evaluated include positive vorticity 
advection, mass divergence and the thermal structure, whereas the absolute vorticity and 
the degree of reduced static stability (d@/dp) of the lower troposphere are also compared. 

This thesis is organized into seven chapters. The second chapter reviews applicable 
literature on explosive cvclogenesis. Chapter Three concentrates on the methods of data 
assimilation and analysis of FGGE data. The separation of explosive and nonexplosive 
cyclones will be described in Chapter Four. The application of the quasi-Lagrangian 
diagnostic technique to this data set will be described in Chapter Five. Chapter Six 
provides a statistical comparison of the two storm groups. The last chapter summarizes 
the results of this research and offers recommendations to improve studies of explosive 


cyclogenesis. 


ll. LITERATURE REVIEW 


The physical processes thought to be responsible for explosive cyclogenesis are la- 
tent or sensible heating, low-level baroclinity, upper-level perturbations and symmetric 
instability. The roles of these events are presented in this chapter by examining current 
literature: 

Sanders and Gyakum (1980) focused attention to the problem of rapid maritime 
cyclone development. Because these low pressure systems form quickly and are fre- 
quently intense thev represent a potentially dangerous environment for mariners. 
Sanders and Gvakum compiled a svnoptic dvnamic climatologv during the period from 
September 1976 to May 1979. Thev concluded that rapid cvclogenesis was primarily a 
maritime event with maximum occurrence in the western oceans. The most frequent 
regions of formation were within or just north of the Gulf Stream and Kuroshio cur- 
rents. In the Aleutian and Icelandic low regions, there was a noticeable lack of explosive 
cvclogenesis formation. One of the primary results of their studv was that these events 
tend to occur in regions with the largest sea-surface temperature gradients. 

The development scenario proposed by Sanders and Gvakum (1980) was one with 
cold polar air moving over the warmer ocean, which results in destabilization of the 
lower atmosphere. This process includes the convection and latent heat release in an 
environment of reduced static stability playing an important role in rapid cyclone de- 
velopment. In addition, the upper-level vorticity advection and low-level baroclinity 
were thousht to be key elements 1m the development precess. 

Uccellini (1986) investigated the role of upper-level baroclinic processes in rapid 
cvclogenesis. An analvsis of the Queen Elizabeth II (QE II) storm of 9-11 September 
1978 revealed the presence of upper level baroclinic processes upstream of the develop- 
ing cvclone. Specifically, a deep frontal zone associated with a short wave trough and 
polar jet were positioned upstream. The baroclinic processes associated with the system 
could have aided in rapid cyclogenesis. The upper-level pattern included a short wave 
trough and a polar jet with maximum winds of 65 m/s and significant vertical and hori- 
zontal wind shear. At the upper level (300 mb), positive vorticity advection and diver- 
gence were present, in addition to an intense frontal zone at the low level (S00 mb). A 


tropopause fold brought upper tropospheric air down to the 700 to $00 mb levels prior 


to the rapid development phase of the cvclone. This study clearly showed that upper- 
level forcing was a key ingredient in the development of the QE II storm. 

A thorough synoptic investigation of frontal wave cyclogenesis in the North Atlantic 
and North Pacific Oceans was conducted by Wash et al. (1988). Two cases of explosive 
cvclogenesis were evaluated using the European Centre for Medium-range Weather 
Forecasts (ECMWF) objective analyses during FGGE. The first case involved the de- 
velopment of a cyclone on a frontal baroclinic zone in the western North Pacific Ocean. 
The incipient cyclone formed due to the surface convergence between an extreme cold 
outbreak and warmer tropical air. Reduced static stability was caused by sensible heat- 
ing of the lower troposphere by a warmer ocean. This resulted in decreasing the braking 
tendency due to vertical motions and upper-level divergence. The second case involved 
the formation of a polar low in the North Atlantic Ocean region. In both cases, strong 
upper-level support was present in the form of an intense jet streak that aided in the 
explosive deepening of the cvclone. The positioning of the divergence region of the 
upper-level jet greatly enhanced the vertical motion and low-level convergence associ- 
ated with the storm development. When the cyclone moved north of the jet in the North 
Pacific Ocean case, the decreased upper-level support resulted in slower development of 
the storm. The results of this study suggested that upper-level forcing was a key element 
in the initiation of explosive oceanic development. This agreed with previous papers by 
Uccellini et al., who concentrated on the role of the upper tropospheric jet and wave 
Structures. However, reduced stability in the lower troposphere certainly aided the 
baroclinic environment necessary to rapidly develop a cyclone. 

In contrast to the Uccellini et al. and Wash et al. analvses. Pagnotti and Bosart 
(1954) emphasized the role of the low-level baroclinic environment and the convective 
bulk heating effects on rapid cyclone development during the Presidents’ Dav storm (18 
February 1979). Significant features of this storm were the boundary laver warming and 
the relatively large moisture flux. The important physical mechanisms in the develop- 
ment of the incipient cyclone were differential heating, moistening and differences in to- 
pographv between the land and sea, in addition to the cold air being dammed against the 
Appalachian Mountains. Rapid deepening occurred when the upper-level southwesterly 
flow was positioned over the cyclone. The convergence of moisture flux in the surface 
to 700 mb layer intensified and expanded eastward. In the coastal boundary laver, 
convergence of the water vapor dominated the horizontal advection of moisture. This 
Was a vital factor in the development of the large amount of precipitation that accom- 


panied this storm. The latent and oceanic sensible heat fluxes demonstrated that the 


C 


boundary laver processes were critical in the development of the storm. The onshore 
wind flow was also a kev element in the resulting frontogenesis and cvclogenesis. This 
provided the additional support for the important ageostrophic circulation as a precur- 
sor for later development. 

The Presidents’ Dav and QE II storms are obviously related to the dvnamic and 
diabatic processes over the entire troposphere. However, the role of upper-level forcing 
can have a variable nature. The position of the upper-level jet streak and the stability 
of the air are critical factors in rapid cvclogenesis. Therefore, the determination of the 
vertical and horizontal structure of the atmosphere 1s important to this evaluation. 

Rogers and Bosart (1986) examined the three-dimensional structure of cvclone de- 


velopment using rawinsonde data. The four stages of the cyclone were defined as: 


e Incipient. The initial formation of a low pressure area as analvzed by the National 
Meteorological Center. This period included the beginning of the cvclone until the 
Start of the rapid deepening rate. 


e Explosive. The period with a deepening rate of at least one mb per hour for 12 to 
24 h. 


e Mature. The period in which the rapid deepening rate ends, and the central pres- 
süre and areal extent of the cvclone's circulation remain quasi-steady. 


e Decaving. During this segment there is a significant increase in the storm's central 
pressure With an associated decrease in intensity. 

This study evaluated the generation of an explosively deepening cyclone that originated 
in a shallow. lower tropospheric system. It subsequently evolved into a deep and intense 
vortex with a strong baroclinic signature. Other contributing factors included jet streak 
dvnamics (Uccellini et al., 1984) and symmetric instability (Emanuel, 1983; Rogers and 
Bosart, 1986). The possibility of nonlinear interactions in the baroclinic instability 
process was also considered. 

One of the methods used to evaluate rapid cyclogenesis was the quasi-Lagrangian 
diagnostic budget technique (Wash et al., 1988). This approach used a spherical coor- 
dinate svstem and transport equations described bv Johnson and Downey (1975a.b). 
The motivation for using a budget analysis is the simplicity of describing the physical 
processes in terms of balancing the mass’circulation inflow and outflow. The inter- 
actions between the budget volume and the environment are analvzed through the 
transport relationships based on the governing equations. This storm-centered budget 
approach provides a framework to composite the explosive and nonexplosive cases in 


this study. 


The summarv of this literature review 1s that the physical processes that accompanv 
explosive cvclogenesis were varied and inconclusive. This thesis will focus on the 
upper-level forcing and the role of low-level baroclinitv in rapid evclogenesis. The spe- 
cific areas of interest are: 

e the nature of the favorable marine environment; 
e the structure of low-level baroclinicity of the incipient cyclone; and 


e the coupling between the upper and lower troposphere. 


The relative magnitudes of these processes will be analyzed in detail to identify the most 
dominant feature in explosive cvclone development. The newest FGGE analvses are 
used to evaluate these elements. The next chapter explains the methods and techniques 


of data assimilation and analvses used in the FGGE processing. 


HI. DATA ASSIMILATION 


The scarcity of data over open-ocean areas limits the study of maritime rapid 
cvclogenesis. The First GARP Global Experiment (FGGE) represented the first time a 
truly global data set has been prepared (Halem et al., 1982). This thesis uses the revised 
FGGE analysis data (Paegle, 1986). The sources for the FGGE data collection included 
surface (sea, land and drifting buov) reports and rawinsondes. Additional data were 
provided from dropsondes, aircraft, pilot balloons and satellite measurements. The re- 
quirement for 500 km horizontal resolution was fulfilled by more than 7000 temperature 
sounding profiles per dav and 6000 cloud drift winds from five geostationary satellites 
(alem era 195». 

The European Centre for Medium-range Weather Forecasts (ECMWTF) prepared 
analvses of the FGGE data set for the entire vear. During the special observing periods 
(SOP s). the archiving was accomplished every 6 h, otherwise it occurred every 2am 
The analvses were performed for 19 levels in the vertical and with a horizontal resolution 
of 1.875° latitude longitude. 

The processing and management of data were separated into three levels of control. 
Level I was comprised of primary data from a variety of instrumented platforms. Level 
II-b represented (where b denotes data collected from the global research data set) the 
transformation of Level I data into basic meteorological parameters (Bjorheim et al., 
1981). The Level II1-b data set contained both basic analvsis parameters in addition to 
derived ones. The original analvses were not initialized and were comprised of 
geopotential height. sea-level pressure and horizontal wind components. The derived 
terms included temperature, relative humidity and vertical velocitv at each vertical level 
(Benetssencet dll 952). 

One of the primary differences between the main and final Level II-b data sets was 
the addition of world-wide synoptic and ship data, which more than doubled the original 
amount of data. The final data set also included additional aircraft reports, data from 
three monsoon experiments, the United States special effort satellite temperatures 
(SATEM s) and cloud-drift winds (SATOB’s), a new set of SATOB winds produced by 
the University of Wisconsin from Japanese HIMAWARI imagery and from the 
METEOSAT for 06 and 18 UTC during SOP’s, and reprocessed DROPSONDE data, 
(Uppala, 1986). 


The initial data assimilation used to produce the Level III-b analvses consists of a 
multivariate optimum interpolation analysis, a nonlinear normal mode initialization, and 
a high resolution model that provided a first-guess forecast. The new ECMWF assim- 
lation of the FGGE Level III-b data consists of the following changes (Shaw et al.. 
1987): 

* diabatic nonlinear normal mode initialization; 

e incorporation of a diurnal cycle and improved humidity analysis; 

e greater vertical resolution in the upper levels; 

e revision of the Optimum Interpolation (OI) statistics package; and 


* arevision of the quality control criteria. 


The present ECMWTF spectral model with rhomboidal truncation and mean 
orographv was used for the FGGE analyses. The revised physical parameterizations 
included a formulation of deep and shallow convection. as well as cloud and radiation 
schemes. Two previous weaknesses in the OI technique were the unrealisticallv low an- 
alvses error in data-rich areas and a high growth rate of forecast error for most regions. 
The solution to this problem has been to assume the revised growth rate is linear, which 
has resulted in agreement with statistical estimates. In addition, the new first-guess es- 
tiniates are lower and sinular to upper-level geopotential values. Another improvement 
in the new FGGE analvses was to increase the vertical resolution from 12 to 19 levéis 
{the horizontal resolution remained at 1.875°). 

The updated quality control methods have introduced checks to eliminate incorrect 
data inputs. Shaw et al. (1987) have described certain observational platforms that have 
consistently provided erroneous data. These platforms have been excluded in the new 
analyses. The ECMWF algorithm for generating the upper atmosphere background 
field for the 10, 20 and 30 mb analvses was improved to include upper-level wave effects. 
This modification was accomplished to include Rossby and Kelvin wave effects and to 
produce more reasonable height fields. 

The comparisons of the initial and final data sets are still being performed. Prelim- 
inary results have shown changes up to ten mb in the Southern Hemisphere and four 
mb in the Northern Hemisphere oceans for a typical mean sea-level pressure pattern. 
The adjustment of wind data in the upper troposphere has been improved while the 
comparison at 400 mb is closer to the initial analysis due to the smaller assumed obser- 
vational errors. Evaluations of ECMWF analyses have shown that greater weighting 


has been given to the more accurate data (e.g., pilot reports and aircraft data) when 


available. In data sparse regions, 1.e., across the open ocean, greater use of SATOB data 
has occurred. 

There have been substantial improvements to the FGGE Level III-b data set 
produced by ECMWF. These improvements and the increased data availability have 
resulted in an appreciably higher qualitv analysis. The locations of rapid cyclogenesis 
and the calculations of storm environment parameters rely heavilv on the accuracy of 
this data set. The progress made in developing the new FGGE data set should improve 
the outcome of this research. The following chapter concentrates on the selection of 


storms and the development of a list of rapid cvclogenesis cases. 
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IV. SCREENING AND SELECTION OF EXPLOSIVE AND 
NONEXPLOSIVE CYCLONES 


The region of cyclone development was based on the studies by Smith (1986) and 
the previous cyclone climatology of the North Atlantic and Pacific Oceans (Sanders and 
Gyakum, 1980). The north-south boundaries (20°N to 60°N) were the same for the 
North Pacific and Atlantic Oceans. The east-west borders in the North Pacific Ocean, 
went from 170°W across the international dateline to 120°E. The east-west boundaries 
in the North Atlantic Ocean stretched from 40°W to 85°W. The period of the evaluation 
was from 17 January to 28 February 1979 as part of the first special observing period 
(SOP). 

The initial screening of the storms was accomplished by Smith (1986) from sea-level 
pressure analvses produced bv the National Meteorological Center. The criterion for 
selection of a cvclone was a closed 1sobar throughout the initial 24 h development pe- 
rod. The requirement for an explosive cvclone was a deepening rate of at least one mb 
per hour for a time period of 12 to 24 h. Cvclones that deepened at a lesser rate were 
classified nonexplosive. The deepening rate was adjusted for latitude bv the correction 
factor sin 0'sin 60%, where O was taken at the point of the most rapid deepening. This 
adjustment factor normalized the pressures for 60^ latitude. In addition, onlv storms 
that deepened over oceanic waters were considered. Storms that regenerated after an 
initial deepening Were not included due to other physical processes being responsible for 
their formation. 

The storm tracks in the North Atlantic and Pacific Oceans are displayed in Figs. 1 
and 2. As previously mentioned, favorable areas of explosive and nonexplosive cyclone 
generation are located near the western ocean boundary currents. 1n the North Pacific 
Ocean, most of the storms occlude and move off to the northeast before reaching the 
international dateline. In the North Atlantic Ocean, the storm tracks are more domi- 
nant off the east coast of the United States, which poses a particularly dangerous threat 
to vessels in the maritime shipping lanes. The most severe storms formed in a 400 to 
500 n mi radius from Cape Hatteras, North Carolina. 

The initial list compiled by Smith (1986) was revised based on the new FGGE sea- 
level pressure analyses. The results of this revision for the North Pacific Ocean are as 


follows: 


® shifted the development positions of the explosive cyclones Pl, P35, لل ل‎ 
the nonexplasive cyclones NPI. NES NPS eee 


e started storms P4, P7, P8, NPI, SP9 and [1نلح‎ 121 
€ started storms NP6 and LAI lao. 
e shifted NP9 and NPI1 to the explosive category, 
e designated P4 and P$ as rapid deepening storms for only 12 h; and 
e dropped NP4, NPS and NP7 from the storm list due to pressure actually increasing 
vice decreasing (discovered in the new FEGGE daa" 
In the North Atlantic Ocean, the following revisions were made: 
e shifted the development locations for storms Al, A3, A4, A6, A7, AS and NA2; 
e started storm A] I? h earlier; 


* dropped NA3 from the storm list due to pressure actually increasing Vice decreasing 
(discovered in the new FGGE data): and 


e noted that A3 development was influenced bv land vice maritime processes. 


Listings of the explosive (A or P) and nonexplosive cyclones CNA or NP ANA 
North Atlantic and Pacific Oceans and their 0, 12 and 24 h positions are provided in 
Tables 1 and 2. The total deepening for the 12 and 24 h periods are listed in Tables 
and 4. In addition, Tables 5 and 6 provide a listing of the changes in the sea-level 
pressure (AP) for each 12 h period. Wherever possible, the center of a closed 1sobar was 
used to begin the initial development stage. In other cases, explosive development oc- 
curred so rapidlv that the 12 h predevelopinent did not show a closed isobar. In those 
cases, rideing in the 500 to 1000 mb thickness pattern and the general structure of the 
surrounding isobars were used to locate the center. 

The results in Tables 1 through 4 indicate that storm development near the western 
ocean boundary currents is quite dominant. The frequency of explosive and nonex- 
plosive cyclone development in this 42-dav sample generally agrees with the previous 
studies (Sanders and Gvakum. 1980) in that the North Pacific Ocean has more storms 
but the North Atlantic Ocean cyclones have a larger deepening rate. Tables | through 
6 list 16 explosive (including the one that formed over land) and the eight nonexplosive 
storms. The cyclones that deepened for only 12 h and the ones that were not formed 
over the ocean environment Were excluded in further statistical analvses. 

The next chapter discusses the development and usage of the mass and circulation 
budgets. These methods will produce the data necessary to determine the most domi- 


nant characteristics in the process of cvclogenesis. 


Table 1. EXPLOSIVE CYCLONE POSITIONS IN THE NORTH PACIFIC 
AND ATLANTIC OCEANS, FROM ECMWF ANALYSES, 1979. * de- 
notes storms with 12 h of deepening, ** denotes initial dev elopment Over 
land and time 1s in UTC. 
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North Pacific Ocean storm tracks and positions for a 24-h interval. 
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Table 2. NONENPLOSIVE CYCLONE POSITIONS IN THE NORTH PACIFIC 
AND ATLANTIC OCEANS, FROM ECMWF ANALYSES, 1979. Time 
E. 
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Table 3. PRESSURE VARIATIONS FOR EXPLOSIVE CYCLONES IN THE 
NORTH PACIFIC AND ATLANTIC OCEANS, FROM ECMWF ANA- 
LYSES, 1979. * denotes storms with 12 h of deepening, ** denotes initial 
development over land and time is in UTC. 
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Table 4. PRESSURE VARIATIONS FOR NONEXPLOSIVE CYCLONES IN 
THE NORTH PACIFIC AND ATLANTIC OCEANS, FROM ECMWF 
ANALYSES, 1979. Time is In A 
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Table 5. PRESSURE DECREASES OF EXPLOSIVE CYCLONES IN THE 
NORTH PACIFIC AND ATLANTIC OCEANS, FROM ECMW F ANA- 
LYSES, 1979. * denotes storms with 12 h of deepening, ** denotes initial 
development over land and time 1s in UTC. 
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Table 6. PRESSURE DECREASES OF NONEXPLOSIVE CYCLONES IN THE 
NORTH PACIFIC AND ATLANTIC OCEANS, FROM ECMWF ANA- 
EYSES, 12708 tame is ange ie 
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V. MASS AND CIRCULATION BUDGET ANALYSIS OF THE CYCLONE 
DEVELOPMENT 


The quasi-Lagrangian diagnostic technique (Wash, 1978; Johnson and Downey, 
1975a,b) is used to evaluate cvclone development, maturation and dissipation for these 
groups of explosive and nonexplosive storms. The primary objective of the storm- 
following method 1s to analyze cyclone development within the westerly wave regimes 
of the midlatitudes. The computations are performed in a spherical coordinate system 
that translates with the storm. This system divides the lateral transport and advection 
associated with the wave translation from the quantities relative to the moving cyclone. 
The proper positioning of the budget volumes is a vital factor in obtaining meaningful 
results from these calculations. 

The volume of the environment is centered on the lowest sea-level pressure from the 
ace analysis. A reference radius 1s constructed from the center of the cvclone to the 
center of the earth. The radius of the volume around the storm is separated into one 
degree latitude increments. This configuration is a cylindrical volume since the radius 
of the earth is much larger than the depth of the atmosphere. 

The advantages of using this method are: 

e focuses in on the kev elements of development for a particular cyclone: 


e evaluates the measure of the parameters responsible for development with respect 
foetne cyclone centers; and 


e quantitatively analyzes the development through the budget and transport re- 
lationship. 

The equations used to describe the mass and vorticity circulation budgets are defined 

in Appendix A. The outputs of these calculations are used to determine the magnitude 


of the following key cyclone development terms: 
e changes in central and areal average pressure; 
e absolute vorticity; 
e vertical velocity; 
e jow-level baroclinicitv; 
e javer-averaged potential temperature; 
e static stability; 


e vorticity advection; 


e eddy mode of vorticity transport; 
e mean mode of vorticity transpor and 


© upper-level divergence. 


The goals of these budgets were to measure the strength of cyclogenesis, and the 
degree of upper- and lower-tropospheric influences. The absolute vorticity, kinematic 
vertical velocity and changes in central and areal pressures are used to evaluate the 
strength of the cvclogenesis. The vorticity advection, eddy and mean modes of the 
vorticity transport, and upper-level divergence are evaluated for upper tropospheric in- 
fluences. The upper-level divergence is measured by integrating the mass transport di- 
rected out of the budget volume in the upper troposphere. Finally, the static stability, 
laver-averaged potential temperature and low-level baroclinity are examined for low- 
tropospheric effects. 

The mean and eddy modes comprise the horizontal! transport of absolute vorticity. 
The lateral transport of absolute vorticity by the eddy mode strongly correlates with the 
Vorticity advection (Calland, 1983). The boundary flux, the mean absolute vorticity and 
the normal component of the mean wind at two adjacent levels Were used to determine 
the laver-average mean mode vorticity. The eddy mode was computed bv subtracting 
the mean mode from the total absolute vorticity transport. The total transport was de- 
fined as the boundary flux of the product of the absolute vorticity and the normal com- 
ponent of the wind (Wash, 1978). 

The kinematic vertical velocity was obtained from the vertically-integrated continu- 
itv equation in pressure coordinates. The divergence was computed by determining the 
Vertical flux of mass at each level (assuming zero at the top and bottom boundaries). 
The lateral flux of mass was calculated from the wind fields and was balanced bv the 
vertical flux to maintain conservation of mass. Unfortunately, errors in wind measure- 
ments resulted in the occurrence of a residual when the horizontal! and vertical fluxes 
were compared. The O’Brien (1970) correction scheme was utilized to distribute these 
residuals. The weighting of this adjustment increased at the upper atmospheric levels 
since the winds were more difficult to measure at higher heights. 

Several methods to measure the strength of the low-level baroclinicity were consid- 
ered. Vertically averaged temperature differences at the storm boundary on a radial 
through the cyclone center were measured. The following different radials were consid- 


ered when evaluating this term: 


e a radial perpendicular to storm motion: 


لم 
tv‏ 


e a radial through the storm center from the warmest area in the warm sector; and 
e a radial through the storm center to the coldest point north or northwest of the 
cyclone. 
The maximum temperature difference from the warmest point along a radial through the 
storm center was chosen as the best technique due to its consistency and reproducibility. 
In general, this radial was the same radial perpendicular to the storm’s motion. The 
largest temperature differences through the cyclone were also found along this radial. 
This temperature difference was a measure of the strength of the low-level frontal zone 
or baroclinity. The static stability was derived from the 1000 to 500 mb temperature 
analyses by the ECMWF. 
Results of the computations of these cyclone development parameters are included 
in the next chapter. The results highlight the major differences in the storm environment 


properties between the explosive and nonexplosive groups. 


25 


VL- RESUET: 


A. INTRODUCTION 

In this chapter a statistical evaluation of the thirteen explosive and eight nonex- 
plosive cyclones was conducted. The two explosive storms that deepened for only 12 h 
and the explosive storm that formed over land were excluded from this statistical study. 
Smith (1986) previously prepared statistics on upper- and low-level cyclone development 
processes. This study uses the new FGGE analyses and focuses on assessing the signif- 
icance of the terms responsible for explosive cyclogenesis. These terms, introduced in 
Chapter Five, were compared for the 24 h cyclogenesis period and at the different time 
intervals (0, 12 and 24 h). The evaluation of these parameters was performed in the 
storm-centered framework using averages from 4^ and 6^ latitude radii budget volumes. 
Hereafter, these volumes will be referred to as 4° / 6° radii volumes. 

Statistical box plots provide a convenient method to display the variations in data 
distribution using quartiles. The change in pressure (AP) shown in Fig. 3 1s used as an 
example. The first two plots (going left to right) are at the 12 h changes in pressure for 
the explosive and nonexplosive storms. The two plots on the right are for the 24 h time 
period. The medians (horizontal line segment within the box) for the 12 h period. are 
11 mb for the explosive cvclones and 5 mb for the nonexplosive ones. The top os 
box for the explosive storms at the 12 h period represents the 75 percentile (15 mb) and 
the bottom is the 25 percentile (8 mb). The interquartile range (Q) is the distance from 
the bottom of the box to the top (Chambers et al., 19835). The value of Q is used to 
define the following terms: 

e adjacent values are within 1.5Q of the quartiles; 


e outside points are between 1.5 to 5Q of the quartiles; and 


e detached points are greater than 3Q away from the quartiles. 


The clear separation of the explosive and nonexplosive 24 h pressure changes in the box 
plots reflects the rapid development definition of the two cvclone groups from the sea- 
level pressure maps. 

The box plots allow a partial assessment of symmetry. If the distribution is sym- 
metric. then the box plot is symmetric about the median. Additionally, if the distance 


to the upper adjacent point from the top of the box is greater than the lower distance 


from the bottom, then the data distribution 1s skewed to the higher pressure, as in the 
case of the explosive storms (for the 12 h mark in Fig. 3). The explosive storms have a 
higher mean (represented bv a circle) and median than the nonexplosive ones, following 
the definition of the groups. The variability of the data for the explosive storms 1s 
greater during both time periods which can be observed in Tables 5 and 6 (Chapter 4). 
This graphical representation provides a method to compare average values and also the 
distribution of the data set. 

The Kolmogorov-Smirnov test for goodness of fit was performed to determine if the 
individual terms are normally distributed (Appendix B contains a more detailed de- 
scription of this process). If the data were normally distributed, then a student t-test and 
90 percent confidence level were calculated to determine the statistical relationships be- 
tween the storm groups. Those storms that met the critical test statistic value for re- 
jecting the hypothesis (that the means were equal) at a significance level of ten percent 
are denoted by a subscript ‘t’ in the following tables. The degree of difference was a 
function of the sample size and the size of the standard deviations (Devore. 1987). The 
standard deviations are included in parentheses next to the mean values in Tables 7 
through I$ for for the upper-level and low-level forcing functions. 

The time tendency and 24 h evolution of upper-level and low-level forcing terms are 
displayed in Figs. 4 through 21. The inner two plots (going left to right) correspond to 
the explosive and nonexplosive storms at 4° latitude radius and the outer plots at 6° ra- 


dius. These figures are for the initial, 12 and 24 h time intervals. 


B. TIME TENDENCY OF LOW-LEVEL CIRCULATION 

The terms that describe the strength of the cvclogenesis process are associated with 
the low-level absolute vorticity and vertical motion. The strength of the low-level cir- 
culation was measured by comparing the absolute vorticity for the 1000 to 850 mb level. 
The initial values for the nonexplosive storms were higher by about ten percent as shown 
in Table 7. Several explosive storms developed so rapidly that no closed isobar are 
present at the initial time. The inclusion of these systems likelv explains the higher initial 
low-level vorticity values in the nonexplosive set and the large vorticitv increases in the 
explosive set. The fact that the initial values were less for the explosive cases indicated 
that the low-level intensity was not the initial dominant factor in the separation of ex- 
plosive and nonexplosive storm development. 

As the storm development progressed, the explosive cvclones exhibited a larger ab- 


solute vorticity increase and final 24 h value averaged over 4° latitude radius. The 24h 


to 
t 


value is the onlv statistically different absolute vorticity between the explosive and non- 
explosive samples, compared to other radii and time periods (Table 7). The higher val- 
ues of absolute vorticity at the initial time for the nonexplosive cyclones suggests that 
Separating the storm groups using these initial values 1s impossible. The greater absolute 
vorticity at the 4° radius for the explosive storms is shown in Figs. 4 through 6 (as time 
increased to the 24 h point). This indicates a tighter pressure gradient and a more or- 
ganized storm in the explosive cases, in agreement with the central pressure statistics. 
The absolute vorticity of the nonexplosive cyclones for the 6° radius volume was 
larger than the explosive storms bv ten percent at the start of the cvclone development. 
During the 24 h development period for the explosive storms, the 6° latitude vorticity 
values do increase more than the nonexplosive storms, but not as dramatically as for the 
4° volume. At the 24 h mark, the average vorticity values are similar for the two groups. 
This indicates that the explosive storms are more compact with an intense inner core, 
but the average vorticity does not differ from the nonexplosive storms using the 6° radius 


volume. 


C. TIME TENDENCY OF VERTICAL MOTION 

The sample statistics for the $50 mb and 700 mb kinematic vertical velocities (units 
of 10-*mb sec) at 4° and 6^ latitude radii are listed in Tables 8 and 9. Figs. 7 through 9 
are the box plots of the vertical velocity for each time interval. The kev results from the 
tables are the statistical separation of vertical motion for the 4° and 6° radu volumes at 
the start of cvclogenesis. This indicates a significant upward vertical motion in the early 
stages of storm development. The maximum vertical velocity occurs at the 12 h point 
for all other radi and time periods. 

Kinematic vertical motion values at 4° radius, for explosive cyclones are greater than 
the nonexplosive ones for all three time periods. The two groups are statistically sepa- 
rable for the O h and 24h periods. After 12 h, the variability of the vertical motion for 
the nonexplosive storms is quite large. This is shown by the sizeable standard deviations 
in Table 8 and the range of values on the box plots in Fig. 8. Consequently, this time 
period did not pass the statistical test. 

The results of the 4° radius vertical motion are also reflected for the 6^ volume. 
Significantly larger kinematic vertical motion values are present for all time periods for 
the explosive storms. The 850 mb vertical motion values are statistically separable for 


all time periods, while the 700 mb vertical motion is separable at the 0 h and 12 h mark. 


The stronger vertical motion values at both the 4° and 6° radii and at most of the 
time periods indicates the explosive set 1s characterized by vigorous upward vertical 
motion that is present at the start of sea-level cvclogenesis. These large vertical motion 
values are consistent with the strong vorticity increases and pressure falls shown earlier. 

In addition to the kinematic vertical velocities, vertical velocities from the ECMWF 
analvses were also studied. Cyclone-centered averages (not shown) from the ECMWF 
vertical motions were considerably more variable in time than the kinematic vertical 
motions. The ECMWF vertical motion results were inconsistent with the behavior of 
the cyclones studied. This suggest inaccuracies in the normal mode initialization of the 


ECMWFE analyses, perhaps caused by a smoothing of the initial vertical velocities. 


D. TIME TENDENCY OF STATIC STABILITY 

The static stability is evaluated in the laver of 1000-500 mb (Tables 10 and 11). The 
results show both explosive and nonexplosive systems develop in a weak static stability 
environment. However, statistical separation between the groups only occurred for the 
6° radius case at the 12 h point. The similarity in the standard deviations and means for 
both storm groups suggest that the low-level environment was basically the same for the 
explosive and nonexplosive cyclones. 

Static stability values (dO, dP) for the explosive storms are 0.2 K/100 mb to 0.5 
K/100 mb less than the nonexplosive storms (Tables 10 and 11). This agrees with the 
expected reduced static stability for the more rapid deepening storms. However, the in- 
itial static stability analyses, hike the initial low-level vorticity, can not be used to sepa- 
rate the nonexplosive from explosive cyclone groups. Although the average static 
stability is similar for the explosive and nonexplosive cyclones, there still may be partic- 
ular quadrants or regions of the cvclone for which there are large static stability differ- 
ences. Quadrant analyses were not conducted in this study but should be investigated 


in the future. 


E. TIME TENDENCY OF LOW-LEVEL BAROCLINITY 

The low-level baroclinic structure is specified from 1000-800 mb and 1000-500 mb 
temperature differences across the frontal zone. Specifically, the temperature values 
through the warmest sector and usually perpendicular to the storm motion are consid- 
ered. The explosive storm temperature differences are slightly larger than their 
counterparts as shown in Tables 10 and 11. However, the differences between the storm 


groups are within a standard deviation, and statistical separation is impossible. 
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The frontal zone strength does not change appreciably during the 24 h development 
period. However, the frontal zone does show a tendency to strengthen for the first 12 
h. The largest gradients for the explosive and nonexplosive cases occurred at the 12 h 
point with the exception of the 6^ latitude radius explosive cvclones at the 1000 to 800 
mb laver (Tables 10 and 11). Overall, the 1000-500 mb values had lower standard devi- 
ations compared to the 1000-800 mb level. This suggests that there is less variability in 
temperature values when including the layer between 800-500 mb. This is consistent 
with the storm-to-storm variations being less when evaluating temperature changes 
through a deeper laver of the troposphere. The results of this term show that the 
strength of the low-level frontal zone does not differ significantly between the explosive 


and nonexplosive groups. 


F. TINE TENDENCY OF UPPER-LEVEL FORCING 

The strength of the upper-level forcing was estimated from the averaged (S00 to 200 
mb) vorticity advection. the eddy and mean modes of vorticity transport. and the 
upper-level divergence. Major differences between rapid and normal cyclogenesis at the 
initial period are found in the upper-tropospheric results. Averaged over 4° latitude ra- 
dius (Table 12), the upper-tropospheric comparison indicates that there 1s stronger pos- 
itive vorticity advection, larger inward eddy mode and outward mean mode vorticity 
transport at the initial and subsequent times. Averaged over 6° latitude (Table 13). the 
positive Vorticity advection is still significantly stronger for the explosive storms. There 
are consistent differences in the 6° radius eddy and mean mode of vorticity transport, 
but not to the degree that allows statistical separation. 

Figs. 10 through 12 show that over a 24h period, the vorticity advection has large 
positive values that vary with time for both storms groups. Specifically, the maximum 
positive vorticity advection was at the 12 h point. The initially strong positive vorticity 
advection became more intense during the first 12 h, presumably through the self- 
development process (Petterssen, 1956). This self-development process is the strength- 
ening of the upper-level wave that results in greater cold advection to the west and 
greater Warm advection to the east. The increase in warm and cold advection act to in- 
tensifv the positive vorticity advection. In the explosive cases, the continuation of 
cvclogenesis after 12 h results in a vertically developed system. The slight decrease in 
magnitude of the values after the 12 h indicates a lessening in the deepening rate of the 


STORMS: 


The eddv and mean mode values are in agreement with the previous discussion of 
the vorticity advection, and peak at the 12 h mark as shown in Figs. 13 and 17 (where 
the inward eddv mode has positive values and the outward mean mode has negative 
values). The eddv mode of vorticity transport 1s strongly related to the average positive 
vorticity advection (PVA) term and its agreement with the PVA term reinforces the im- 
portance of strong upper-level influences on explosive cyclogenesis. The mean mode of 
vorticity transport isolates the role of the upper-level divergence in the cvclone’s vorticity 
budget. The stronger negative mean mode results documents stronger divergence aloft 
for the explosive cvclogenesis set 1n agreement with the vertical motion results presented 
above. 

The upper-level divergence, as estimated by the outward mass transport values, has 
simular fluctuations with significantly larger values for the explosive cases (up to a factor 
of four) compared to the nonexplosive ones (Figs. 19 to 21). The divergence term 
showed the greatest variability of all the upper-level properties, particularly for the ex- 
plosive svstems. The primary reason for this variability was probably the difficulty in 
analvzing the tropospheric winds. 

The upper-level forcing terms for explosive and nonexplosive storms contained sta- 
tistical separation in over half of the categories. The largest separation is at the storm 
onset where the explosive storms had three times the divergence of the nonexplosive 
ones at the 4° radius. The differences between the explosive and nonexplosive cases 
gradually decreased to 50 percent as the deepening progressed. The best separation be- 
tween the explosive and nonexplosive sets of storms for the upper-troposphere parame- 
ters was at the J^ radius. This suggests the important upper-tropospheric features are 


near the surface cyclone center and not spread over a large area. 


G. TWENTY-FOUR HOUR AVERAGE STORM COMPARISON 

The 24-h averaged values used to describe cyclone development are listed in Tables 
14 through 18. The overall strength of the low-level baroclinic, low-level circulation and 
upper-level forcing for the explosive and nonexplosive cyclones can be compared in these 
tables. Table 14 shows the sea-level pressure for both types of storms. The sea-level 
pressure average for 24-h average was 18 mb greater for the explosive cases compared 
to the nonexplosive cases. 

The low-level baroclinic structure of the explosive and nonexplosive storms were 
very simular with the exception of the former having a higher area-averaged potential 


temperature by 4° (for the 4° and 6° latitude radius). The static stability values were less 


(Table 15) for the explosive storms by a margin of approximately 0.5 degrees, which was 
not statisticallv significant. 

]he computed kinematic vertical velocities in Table 16 were consistent with the ob- 
served upper-level parameters. The explosive cyclones have consistently stronger up- 
ward motion at the 700 and 850 mb levels. The absolute vorticity averaged for 24 h 
showed higher values at the 4° radius compared to the 6° radius position, even though 
the initial values were higher for the nonexplosive cyclones (Table 17). However, as time 
progressed, the explosive storms gained intensity and eventually showed a sizeable in- 
crease over the nonexplosive storms. 

The upper-level diagnostic terms averaged over a 24 h period (Table 18), indicate 
that there was a significant amount of of upper-level forcing in the explosive cyclogenesis 
development The vorticity advection was positive in all cases with the explosive storms 
having nearly double the values of the nonexplosive cases (Table 18). The upper-level 
divergence was approximately two to three times larger for the explosive storms (how- 


ever, the amount varied with time). 


H. SUMMARY 

The kev result of this chapter is that the kinematic vertical velocity and the upper- 
level forcing mechanisms are statistically separable. The large values for the upper-level 
processes suggests that the upper-tropospheric wave influence 1s most likelv to produce 
the stronger vertical motions that accompany explosive cvclogenesis. In contrast, the 
low-level baroclinity and static stability did not show a statistical separation. which 1n- 
dicates that the explosive and nonexplosive storms share a similar low-level frontal 
strength in a weak static stability environment. The statistical separation of the ex- 
plosive and nonexplosive cvclones is stronger than in the early results of Smith (1986). 
The reasons for their improvement likely stem from the revised FGGE analyses as well 
as a better delineation of the storm tracks and intensity. 

The evolution of the upper-level developmental parameters is characterized bv 
maximum values occurring at 12 h in the middle of the explosive deepening period. At 
this decisive point, the storm either continued to develop as in the explosive cases or 
Starts to weaken. However, the low-level circulation continues to increase to 24 h and 
does not reach a peak value after 12 h. In addition, the static stability continues to de- 


crease but dt a lesser rate during the lachi? am 


30 


In the next chapter these events are explored and evaluated. The conclusions are 
discussed and recommendations were offered in an effort to gain insight on the impor- 


tant mechanisms of explosive cvclogenesis. 
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Table 7. TWENTY-FOUR HOUR TIME EVOLUTION OF LOW-LEVEL AB- 
SOLUTE VORTICITY FOR BOTH STORM GROUPS (4? AND 6? 
RADI). Standard deviations are listed in parentheses and ‘t’ subscript 
indicates explosive and nonexplosive sets are statistically different at a 
1095 significance level. 


SIZE RADIUS 4? RADIUS 6º 


ABSOLUTE EX- EX- 
VORTICITY(I0E-6/S) | PLOSIVE { NON-EXP | prosive | NONEAF 


ENE 










Table 8. TWENTY-FOUR HOUR EVOLUTION OF KINEMATIC VERTICAL 
MOTION FOR BOTH STORM GROUPS (4° RADIUS). Standard de- 
viations are listed in parentheses and ‘t’ subscript indicates explosive and 
nonexplosive sets are statistically different at a 10% significance level. 















iss qo) a 
SEX 


Table 9. TWENTY-FOUR HOUR EVOLUTION OF KINEMATIC VERTICAL 
NIOTION FOR BOTH STORM GROUPS (6° RADIUS). Standard de- 
viations are listed in parentheses and ‘t’ subscript indicates explosive and 
nonexplosive sets are statistically different at a 10% significance level. 


ONEX- 
DESCRIPTION EXPLOSIVE ee E 


T -156 97) 
=12n “165, (109) 











= 
| | || 
i 、 
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Table 10. TWENTY-FOUR HOUR TIME EVOLUTION OF THERMAL 
STRUCTURE FOR BOTH STORM GROUPS (4? RADIUS). Standard 
deviations are listed in parentheses and 't' subscript indicates explosive 
and nonexplosive sets are statistically different at a 10% significance 
level: 
















DESCRIPTION EXPLOSIVE DUE 
STATIC STABILITY (d0/dP) MEE e | 
I ( * K/100 MB) (° K/100 MB) 


“170 36 —[ 529 (49 
“4.45 (083) | 4.91 (0.38) 
T=24h —4.41 (0.66) —4.61 (0.88) 










STRENGTH 
is 56 
13 (5:31) 





18 (4.39) 16 (4.57) 
15 ( 


5099, 1555 
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Table 11. TWENTY-FOUR HOUR TIME EVOLUTION OF THERMAL 
STRUCTURE FOR BOTH STORM GROUPS (6° RADIUS). Standard 
deviations are listed in parentheses and ‘t’ subscript indicates explosive 
and nonexplosive sets are statistically different at a 10% significance 
level. 


DESCRIPTION EXPLOSIVE a 
STATIC STABILITY (d0/dP) 
1000-500 MB ( ^ K/100 MB) (* K'100 MB) 


ESTE nn 





Urð 
V 


Table 12. TWENTY-FOUR HOUR EVOLUTION OF UPPER-LEVEL PROP- 
ERTIES AVERAGED FROM 500 TO 200 MB (4° RADIUS). Standard 
deviations are listed in parentheses and ‘t’ subscript indicates explosive 
and nonexplosive sets are statistically different at a 10% significance 
level 


> NONEX- 
DESCRIPTION EXPLOSIVE PLOSIVE 


MEAN MODE OF VORTICITY TSPT 
EXE 
Hein RO 

UPPER-LEVEL DIVERGENCE 
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Table 13. TWENTY-FOUR HOUR EVOLUTION OF UPPER-LEVEL PROP- 
ERTIES AVERAGED FROM 500 TO 200 MB (6° RADIUS). Standard 
deviations are listed in parentheses and t subscript indicates explosive 
and nonexplosive sets are statistically different at a 10% significance 
level 


84 (55) 

MEAN MODE OF VORTICITY TSPT 
IES 
372 (666 
825, (685) 
30 (380) 
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Table 14. TWENTY-FOUR HOUR AVERAGE OF SEA-LEVEL PRESSURE 
FOR BOTH STORM GROUPS (4° AND 6° RADII). 


: ; NONEX- 
DESCRIPTION EXPLOSIVE PLOSIVE 
STORM TYPE PRESSURE PRESSURE 
Initial Surface Pressure (MB) 1006 1006 


24 H Surface Pressure Change (MB) 







Table 15. TWENTY-FOUR HOUR AVERAGE THERMAL STRUCTURE FOR 
BOTH STORM GROUPS (4° AND 6° RADII). 


SIZE RADIUS 4° RADIUS 6° 
STATIC STABILITY me EX- l 
(d0/dP)( © K/100MB) Mns ليد‎ | PLOSIVE OD 
cues» xp — TL d5 
— —Á ZONE 
THERMAL STRENGTH 
(AT) 











Table 16. TWENTY-FOUR HOUR AVERAGE KINEMATIC VERTICAL MO- 
TION FOR BOTH STORM GROUPS (+° AND 6° RADII). 


RADIUS 4º RADIUS 6º 


EX- T EX- 
DESCRIPTION PLOSIVE NON-EXP PLOSIVE NON-EXP 


OMEGA sd OMEGA . . | (10E- MB Sy (IOE-JAMB Sy (IOE-AMB S. — 4MB S 
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Table 17. TWENTY-FOUR HOUR AVERAGE OF LOW-LEVEL ABSOLUTE 
VORTICITY FOR BOTH STORM GROUPS (4° AND 6° RADII). 


SIZE RADIUS 4° RADIUS 6° 


P i EX- 
DESCRIPTION E OAE bí vr | NON-EXP 
" "M 6/SEC)| (10E-6/SEC)| (10E-6/SEC)| (10E-6/SEC) 







Table 18. TWENTY-FOUR HOUR AVERAGE UPPER-LEVEL FORCING 
PROPERTIES FOR BOTH STORM GROUPS (500-200 MB AT 4? 
AND 6° RADII). 


SIZE RADIUS 4° RADIUS 6 ° 


EX- EX- 
DESCRIPTION PLOSIVE NON-EXP PLOSIVE NON-EXP 


VORTICITY 
ADVECTION 
(10E-11/SEC2) 


EDDY MODE 
(10E-11/SEC2) 


MEAN MODE 
(10E-11/SEC2) 


UPPER-LEVEL 
DIVERGENCE(I0E- 10 
GM/SEC) 
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Figure 3. Pressure changes after 12 (inner plot) and 24-h for explosive and nonex- 


plosive storms in box plot format. 
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2 3 
TYPE OF STORM (EXPLOSIVE VS NONEXPLOSIVE) 


QUARTILES 1-EXP AT 4 DEC 


MEANS E 
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Figure 4. Average absolute vorticity (1000-850 mb) box plots for initial explosive 


and nonexplosive storms. 
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Figure 5. Average absolute vorticity (1000-850 mb) box plots for the 12-h mark. 
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Figure 6. Average absolute vorticity (1000-850 mb) box plots for the 24-h mark. 
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Figure 7. Kinematic vertical velocity (850 mb) box plots for initial explosive and 


nonexplosive storms. 


12 HOUR POSITION AT 4 AND 6 DEGREES RADII 


—100 0 


—200 


حم 
O‏ 
ở‏ 
Ð‏ 
= 
Ý+‏ 
| 
uJ‏ 
© 
= 
5 
9 
لا 
> 
= 
< 
© 
br‏ 
> 


— 300 


2 ð 
TYPE OF STORM (EXPLOSIVE VS NONEXPLOSIVE) 


Ð †—EXP AT 4 DEC 
2-NONEXP AT 4 DEG 
ADJACENT VALUES 3-EXP AT 8 DEG 


OUTSIDE POINTS 4—NONEXP AT 6 DEG 
DETACHED POINTS 





Figure 8. Kinematic vertical velocity (850 mb) box plots for the 12-h mark. 
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Figure 9. Kinematic vertical velocity (850 mb) box plots for the 24-h mark. 
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Figure 10. Upper-level vorticity advection (500-200 mb) box plots for initial ex- 


plosive and nonexplosive storms. 
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Figure 11. Upper-level vorticity advection (500-200 mb) box plots for the 12-h 


mark. 
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Figure 12. Upper-level vorticity advection (500-200 mb) box plots for the 24-h 
mark, 
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Figure 13. Averaged eddy mode of vorticity transport (500-200 mb) box plots for 


initial explosive and nonexplosive storms. 
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Figure 14. Averaged eddy mode of vorticity transport (500-200 mb) box plots for the 
]2-h mark. 
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Figure 15. Averaged eddy mode of vorticity transport (500-200 mb) box plots for the 
24-h mark. 
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Figure 16. Averaged mean mode of vorticity transport (500-200 mb) box plots for 


initial explosive and nonexplosive storms. 
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Figure 17. Averaged mean mode of vorticity transport (500-200 mb) box plots for 


the 12-h mark. 
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Figure 18. Averaged mean mode of vorticity transport (500-200 mb) box plots for 
the 24-h mark. 


55 


O HOUR POSITION AT 4 AND 6 DEGREES 


500 1000 1500 2000 


0 


on™, 
Oo 
ليا‎ 
un 
> 
= 
U 
O 
| 
لنا‎ 
© 
Vo 
LJ 
© 
zZ 
لنا‎ 
to 
BE 
Lu 
> 
e 
: 
| 
K 
LU 
C 
að 
= 


—500 


2 
TYPE OF STORM (EXPLOSIVE VS NONEXPLOSIVE) 


QUARTILES 1-EXP AT 4 DEC 
MEANS 2-NONEXP AT 4 DEG 
ADJACENT VALUES — 3 gxp AT B DEG 
OUTSIDE POINTS A-NONEXP AT 6 DEG 
DETACHED POINTS 





Figure 19. —. Upper-level divergence, as estimated by the outward mass flux, (500-200 


mb) box plots for both storm groups. 
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Figure 20. Upper-level divergence, as estimated by the outward mass flux, (500-200 


mb) box plots for the 12-h mark. 
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Figure 21. Upper-level divergence, as estimated by the outward mass flux, (500-200 


mb) box plots for the 24-h mark. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


This study of explosive cvclogenesis using the new FGGE data resulted in several 
conclusions, as follows: 


e the upper-level parameters (statistical separation in over half of the terms) were 
dominant in explosive cvclogenesis, especially during the transition between the 
incipient and explosive phases. 


e the low-level vorticity for the explosive cases displayed a larger growth in magni- 
tude for a 24 h period, in agreement with the sea-level pressure analvses. Even 
though the nonexplosive cyclones had slightly higher initial values (of low-level 
absolute vorticity), eventually the explosive storms developed the larger vorticities. 


e the low-level baroclinity (as measured from a cross-front temperature difference) 
had a maximum value at the 12-h mark with slightly larger values for the explosive 
storms. However, there is no statistical separation in comparing the two storm 
groups. 


e the explosive and nonexplosive storms form in an environment of weak static sta- 
bility. 


e kinematic vertical velocities were more consistent with upper-level physical mech- 
anisms and had less variability than the vertical velocities provided by the ECM WF 
analvses. 


® positioning of the cyclones was critical for obtaining reasonable values for the mass 
and circulation budgets. 


e the new FGGE analvses changed the classification of several explosive and nonex- 
plosive systems and enhanced the overall data set. 


The interpretation of the results listed in Chapter Six stem from the vigorous 
cvclogenesis that was occurring from the onset of the storm to the 12-h point. During 
this period, the cvclone underwent a self development phase in which the strengthening 
of the upper-level wave resulted in greater cold advection to the west and greater warm 
advection to the east. The low-level absolute vorticity continued to grow until the 
pressure of the cyclone started to rise, which marked end of the cyclogenesis process. 

The similarities and differences between the explosive and nonexplosive cyclones 
remained consistent throughout the rapid development in most cases. The maximum 
contribution for the upper and lower processes normally occurred at the 12 h mark. For 
the nonexplosive cases, there was usually a significant decrease after this time. However, 
the nonexplosive storms did experience higher values for absolute vorticity values than 


the explosive storms at the initial cyclogenesis. Therefore statistical separation was not 
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possible for low-level absolute vorticity in the early stages of cyclogenesis. The explosive 
storms have slightly weaker static stability compared to the nonexplosive ones. 
To improve present research efforts the following recommendations are provided: 


e more cases of upper-level forcing are necessary to improve the statistical data base 
for comparison. 


e larger sample sizes of cyclones are needed to provide a better Kolmogorov-Smirnov 
test for goodness of fit for a normal distribution. 


e storms should be grouped for onlv 12 or 2d h of deepening, rather than combined, 
to maintain consistency of the data set. Given the inherent inaccuracies in central 
pressure values, a change in pressure divided by the change in time (AP/AT) has 
more signal (less contamination from ‘noise’) if the total time change (AT) is 24h 
vee l2 


e storms that formed over land and then moved seaward should be excluded from the 
list. Ihis would preserve the purity of the data set as a study of maritime 
cvclogenesis. 


e more detailed studv of upper-tropospheric potential vorticitv is necessary to ob- 
serve upper-level air flow patterns associated with explosive cvclogenesis. 


e different quadrants of the cvclone should be investigated for the largest static sta- 
bility differences. 


e more analyses of the moisture variables would be an important addition O 
processes being evaluated. 

The intensity of the phvsical mechanisms responsible for explosive cyclogenesis are 
varied and, even with modern techniques, difficult to consistently forecast. Although the 
improvement of model skill has occurred due to better analyses, increased resolution and 
treatment of boundary layer fluxes, initial data input is still the most limiting factor in 
predictabilitv (Sanders, 1987). Present methods to obtain the necessarv data to forecast 
these mesoscale events are inadequate. Therefore, the mariner must be knowledgeable 
of the indications for explosive cyclone events while meteorologists continue to pursue 


this challenge. 
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APPENDIX A. BASIC EQUATIONS FOR QUASI-LAGRANGIAN 
BUDGET STUDIES 


l. Definition of Budget Integrated Property (F) 


F= | pJ,fr° sin BaV, (1) 
V. 


where: 
F - is the total property. 
V,- volume of cylinder represented in coordinate system 5 ,o,and f. 
p - densitv. 
| J, | - Jacobian transformation from 5 to z coordinate svstem. 
f - any specific property. 
r - radius in the horizontal. 
B- angle between the local vertical and the edge of the storm measured at the center 
Othe earth. 


2. Mass Budget Equation 


011 > à 
D c LT + FT (2) 
Ahere: 
LT - is the lateral transport of the mass of F. 


VT - 1s the vertical transport of the mass of F. 
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3. Lateral Transport (in detail) 





4, (22 
LT == | | py, (u — w)gfr sin Bdady la, (3) 
na~ O 
My ene: 
7,- 18 the top coordinate. 
s - 15 the bottom coordinate. 
u - relative wind vector in the radial direction. 
Ww - storm motion vector in the radial direction. 
4. Vertical Transport (in detail) 
. dn po. 
VT اح‎ | | 00 men Vx“ sn fdodf | (4) 
o Jo dt dt 1 
Where: 
dy : Å , ; 
uc the change 1n the » coordinate system with respect to time. 
di]; : : 
cmm the change in the top and bottom boundary of the n coordinate system 
£ 


with respect to time. 


5. Vorticity Budget Equation 
The vorticity budget equation shows the contributions to vorticity changes from 
the lateral and vertical fluxes (transports) as sources’sinks. The following equation de- 


scribes these terms: 
C 
S FF LIFT) (5) 


where: 
C, - 1s the absolute circulation. 
6, - 1s the absolute vorticity. 


LT(£,) - 1s the lateral transport of absolute vorticitv that includes the mean and eddy 
modes, as well as the horizontal divergence and advection. 


DVT(£,) - is the vertical transport of absolute vorticity that includes vertical diver- 
gence and advection. 


S(¢,) - is the sink of absolute vorticity through the divergence, tilting and frictional 
dissipation terms. 


The source or sink terms originate from the generation or dissipation of vorticity 
within the budget volume, and in isobaric coordinates represent the divergence, tilting 
and frictional terms. To gain a better physical understanding of the vorticity budget 


equation, the absolute circulation is written in terms of Stokes’ line integral theorem: 


Ca= | t «atm | [v taa (6) 

Ë A 

The total flux u£, can be divided into the divergent and advective components as follows: 
Vel uHt(Veu)tueVe, (7) 


where the terms on the right side of the equation are defined as: 
C(V * u) - vorticity divergence. 


u * VC, - vorticity advection component. 
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APPENDIN B. STATISTICAL ANALYSIS FOR NORMAL 
DISTRIBUTION 


The statistical analvsis for normal distribution on the FGGE cyclone diagnostic 
terms was accomplished using the Kolmogorov-Smirnov test for goodness of fit 
(Lilliefors, 1967). For a sample size », an empirical distribution function E (x) is defined 


as follows: 
F„(x) = kịn (8) 


where: 
P(X) 1s the €mpirical distribution funcion or AN 
k - is the number of observations not greater than x. 


x) 19096 - denotes tlie Sample Values. 


In this test. the null hvpothesis was that the sample had been drawn from a normal 
distribution F,(x), which should be close to F(x). The Kolmogorov-Smirnov test statistic 
D was based on the maximum difference between the observed and hvypothesized cu- 


mulative distribution functions as follows: 
D = max |F,(x) — F(x) (9) 


Comparison of the critical values of D for a two-tailed test was conducted under the 
null hypothesis at a significance level equal to 0.10. The null hypothesis was rejected if 
D exceeded the critical value. This test assumed that the observed distribution was un- 
known and continuous, and that x, through x, did not directly depend on the choice of 
the distribution (Conover, 1980); 

Examples of the graphical representations of the data distributions are shown in 
Figs. 22 through 27. The data for the upper-level divergence (at the 4° latitude radius 
and T=0 h) of an explosive storm group are displayed in the first three figures , and for 
the nonexplosive case in the latter three diagrams. The histograms in Figs. 22 and 25 
illustrate the relative frequency distribution for the 13 explosive and eight nonexplosive 
storms with a normal distribution curve for comparison. Figs. 23 and 26 show the 
normal cumulative distribution function, where the smooth curve represents the contin- 


uous area under the normal distribution curve. The Kolmogorov-Smirnov boundaries 
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are the two diametrically opposed dotted lines and the middle line represents the ob- 
served data. The probability plots (Figs. 23 and 27) demonstrate the same normal dis- 
tribution using a vertical percentile scale. This scale has been modified in the vertical 
to allow an easier visual comparison of the normal distribution (the straight line) to the 


observed data. 
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O HOUR POSITION AT 4 DEGREES RADIUS 
NORMAL DENSITY FUNCTION, N=13 
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Figure 22. Histogram distribution for upper-level divergence for an explosive storm 


group. 
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O HOUR POSITION AT 4 DEGREES RADIUS 
NORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=13 
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Figure 23. Normal cumulative distribution curve for an explosive storm group. 
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O HOUR POSITION AT 4 DEGREES RADIUS 
NORMAL PROBABILITY PLOT, N=13 
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Figure 24. Probability plot for an explosive storm group. 
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O HOUR POSITION AT 4 DEGREES RADIUS 
NORMAL DENSITY FUNCTION, N=8 
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Figure 25. Histogram distribution for upper-level divergence for a nonexplosive 


storm group. 
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O HOUR POSITION AT 4 DEGREES RADIUS 
NORMAL CUMULATIVE DISTRIBUTION FUNCTION, N=8 
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Figure 26. Cumulative distribution curve for a nonexplosive storm group. 


70 


O HOUR POSITION AT 4 DEGREES RADIUS 
NORMAL PROBABILITY PLOT, N=8 
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Figure 27. Probability plot for a nonexplosive storm group. 
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